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Abstract 

It is known that a Green's function-type condition may be used to derive rates 
for approximation by radial basis functions (RBFs). In this paper, we introduce a 
method for obtaining rates for approximation by functions which can be convolved 
with a finite Borel measure to form a Green's function. Following a description of 
the method, rates will be found for two classes of RBFs. Specifically, rates will be 
found for the Sobolev splines, which are Green's functions, and the perturbation 
technique will then be employed to determine rates for approximation by Wend- 
land functions. 

1 Introduction 

Radial basis function (RBF) approximation is primarily used for constructing approx- 
imants to functions that are only known at discrete sets of points. Some advantages of 
this theory are that RBF approximation methods theoretically work well in arbitrarily 
high dimensional spaces, where other methods break down, and ease of implementa- 
tion. Even though the theory tells us that there is no dimensional bound on the appli- 
cability of RBF techniques, in practice, interpolating large sets of data coming from 
high dimensional spaces can be computationally expensive. Fortunately, efficient al- 
gorithms HE] ED have been developed to compensate for such difficulties and allow 
for wider application of RBF methods. As a result of their favorable properties and 
efficient implementation schemes, RBF techniques are being used to solve a variety 
of applied problems, some examples being problems in data mining (4j [SJ, statistical 
learning theory [16], and numerical partial differential equations. In particular, Flyer 
and Wright have shown that RBFs perform well in problems from mathematical geo- 
sciences, l9l [l0l . 

The initial approach to scattered data approximation on W 1 concerned the station- 
ary setting, with RBFs being scaled to be proportional to the fill distance. In this paper, 
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we will study the alternative non-stationary setting, which is relatively new and still 
contains many open problems. For such an approximation, one approach to proving 
error estimates is to exploit a Green's function representation. For example, this tech- 
nique was used by Hangelbroek for analyzing approximation by thin-plate splines on 
the unit disk, lfT4l|T5l . More recently, DeVore and Ron derived rates for the approxi- 
mation of LP functions by RBF spaces Sx(<&), 0- Here, Sx($) denotes a collection 
of linear combinations of translates of an RBF <t>, where the translations come from 
a scattered discrete set of points X C M. d . In order to derive these rates, the functions 
<t> were required to satisfy a Green's function-type condition. The focus of this paper 
will be to extend a result of DeVore and Ron to include certain non-Green's functions. 
Specifically, the Wendland functions, which have some nice properties, such as com- 
pact support, and are hence popular for applications (0 [T3] , do not fit the framework 
of 0. Therefore, this paper will address a generalization of a result from [5 1 that will 
allow us to derive approximation rates for the Wendland functions. Additionally, rates 
for approximation by Sobolev splines will be determined. 

To begin we will prove approximation rates for functions that can be convolved with 
a finite Borel measure to form a Green's function. The motivation for this approach is 
the success of Mhaskar, Narcowich, Prestin and Ward in deriving rates for perturba- 
tions of Green's functions on the sphere, [ 18 1. In fact, they were able to find rates for 
the Wendland functions restricted to the sphere. After establishing this general result, 
we will provide some examples. To be precise, polynomial reproducing functionals 
will be used to show that the Sobolev splines and Wendland functions can be well ap- 
proximated by local translates of themselves, and an error bound for approximation by 
some particular Sobolev splines and Wendland functions will be derived. 

1.1 Notation and Definitions 

The Fourier transform and Laplace transform are fundamental tools for proving many 
results in approximation theory. We will use them to characterize the RBFs under 
consideration and to prove several results. Throughout this paper, we will use the 
following conventions. Given a function / : M. — > R in L 1 , its Fourier transform / will 
be defined by 



Given a function / : [0, °°) — > R that grows no faster than e c " for some a > 0, its Laplace 
transform will be defined by 



Let £ be a subset of W 1 , then we denote by C k (E) the collection of real valued 
functions defined on E that have continuous partial derivatives up to order k. The set of 
functions in C*(R) that converge to at infinity will be denoted by Cq, and we represent 
the compactly supported elements of C k (M. d ) by C*. 

The Schwartz class 5? of functions on Mr is defined as follows. A function / : 
W 1 — s- R is said to be of Schwartz class if for all multi-indices a and /3, there exists a 
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constant C a g > such that 



x a Dpf(x) 



<C, 



forallxeM c/ . 

We shall use the standard definition for the LP spaces. A Lebesgue measurable 
function / : W 1 -> K is in L p {R d ) for 1 < p < oo if 

1/P 



ljf(x)\ p dx) <oo 



and a Lebesgue measurable function / is said to be in L°°(R d ) if 

= ess sup |/(x)| < oo. 



Additionally a measurable function / is said to be locally integrable, denoted by L| oc (M. d ) 
if 

j E \f( X )\dx<°o 

for each bounded, measurable set E C W 1 . 

The space of finite Borel measures on W l will be denoted by M(M c/ ). This space is 
equipped with a norm defined by 

HmINImIO^), 

where \ is the total variation of ji. Given a measure ji e M(M rf ) and an L'^ function 
/, their convolution is an LP function, and it satisfies the following generalization of 
Young's inequality, [ 12 Proposition 8.49]: 

H/*Mllp<ll/UlMl|. (1) 

In deriving error estimates for the Wendland functions, we will need to consider a 
decomposition of measures. Specifically, any measure ji e M(M'') can be written as 
fX = Ha + fl s + fid' where fi a is absolutely continuous with respect to Lebesgue measure, 
fid is a countable linear combination of Dirac measures, and fl s = ji — ji a — jj.^ is the 
singular continuous part of ji. 

Our primary focus will be approximating functions that lie in subspaces of LP 
spaces. The spaces that we will be mainly interested in are the Bessel-potential spaces 
L ' P (M. ), which coincide with the standard Sobolev spaces W k ' p {R d ) when k is a pos- 
itive integer and 1 < p < oo, cf. |fl9l Section 5.3]. The Bessel potential spaces are 
defined by 

L k - p = {f:f=(l + \\-\\ 2 2 )- k / 2 g lg eL p (R d )} 
for 1 < p <°°, and they are equipped with the norm 

H/lb,, = lklL- 
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We will also be working with smoothness spaces associated with linear operators. 
If T : C*(M c/ ) -> C c (R d ) is linear, then we define a semi-norm and norm on C"(R rf ) by 

\f\w(LI'(R d ),T) = \\ T f\\p 
\\f\\w{LP(R d )J) = \\f\\p + \f\w{U'(R d ).T) ■ 

The completion of the space C™(W l ), with respect to the above norm, will be denoted 
by W(L p (M. d ),T). In particular, notice that for any positive integer n and 1 < p < °°, 
we have 

W(L p (K. d ),(l - A)") =L 2n ' p (M. d ) 

where A denotes the d— dimensional Laplacian. 

The approximants used in this paper will be finite linear combinations of translates 
of an RBF <t>, and the translations will come from a countable set X C M. d . The error 
of this approximation, which is measured in a Sobolev-type norm, depends on both the 
function <t> and the set X. Therefore, given an RBF <t> and a set X, we define the RBF 
approximation space Sx (&) by 

S x (4>) = | £ a ^(- - £) : Y C X,#T < °° j r\L l (R d ). 

Note that the functions <t> and sets X will be required to satisfy certain conditions so 
that we may prove results about rates of approximation. The error bounds will be stated 
in terms of the fill distance 

h x = sup inf ||x-€|| 2 , 

xeM d ^ eX 

which measures how far a point in W 1 can be from X, and we additionally require the 
separation radius 

qx = r inf 
2 

to be bounded in order to prevent accumulation points in X. We will therefore be 
working with sets for which the mesh ratio px '■= h-x/qx is bounded by a constant; 
such sets are called quasi-uniform. 

1.2 Radial Basis Functions 

The two classes of RBFs that we will be interested in are the Sobolev splines and the 
Wendland functions. The Sobolev Splines Gy : W 1 — > R form a class of radial basis 
functions most easily defined in terms of their Fourier transforms: 

(5 r (fl)):=(l + ||fl)||lr^ 

for 7 > 0. Notice that Gy is the Green's function for the (pseudo-)differential operator 
(1 — h.yl 2 . These functions have been studied extensively in (HE). I n these papers, it 
is shown that Gy is an analytic function except at 0, and several useful asymptotic ap- 
proximations are determined. Those necessary for our results are listed in the following 
proposition. 
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Table 1 : Examples of Wendland functions 



Function 




Smoothness 


*i,o(*) = ( 


\-r) + 


C u 


= ( 


l-r) 3 + (3r+l) 


c 2 


*1,2W = ( 


1 -r)^(8r 2 +5r+l) 


c 4 


*3,oW = ( 




c u 


*3,lW = ( 


l-r)«(4r+l) 


c 2 


®3,2{x) = ( 


l-r)<j_(35r 2 + 18r + 3) 


c 4 


$3,3 (*) = ( 


l-r)^(32r 3 +25r 2 + 8r+l) 


c 6 



Proposition 1.1. (0 p.253]) For any multi-index a, there are positive constants C 
such that 

for \a\> y—d and \a\ odd 
for \a\ > y—d and \a\ even 
for \ct\< y—d and \a\ odd 
for \a\ < y—d and \a\ even 

Next, the Wendland functions compose a class of compactly supported RBFs that 
are radially defined as piecewise polynomials, and some examples are provided in Ta- 
ble Q] Wendland's book [20 1 (particularly chapter 9) provides a detailed analysis of 
these functions and some of their approximation properties; however, an exact form 
for the Fourier transforms of these functions has not previously been found. Since we 
will need this information for our approximation analysis, we begin by computing the 
Fourier transforms of these functions. The Wendland functions Qdk are determined 
by a dimension parameter d and a smoothness parameter k, and they lie in C 2k (W l ). 
In Table Q] and in what follows, r will be used to denote ||jc|| 2 , and = will be used to 
indicate equality up to some positive constant factor. 

We will now derive an explicit form of the Fourier transform of <I>^ j in the case d 
is odd. Using the notation of ll20l Section 10.5], let d = 2n + 1 and m^n + k. Then by 
ll20l Lemma 6.19] and the definition of O^ i, we have 

&d,k(x)=B m f m (r)r- 3m - 2 (2) 
where B m is a positive constant and the Laplace transform of f m satisfies 

■^fm{ r ) = r m+l(l +r 2)m+l ■ 

In order to find the inverse Laplace transform of the above expression, we will make 
use of partial fractions. First, note that there exist constants a ; -, j3,, and yj such that 

I m a- m B m v- 

(1 + s 2 )m +i L s j+i f L (j + t L (j _ iy+i > > 



Cy 



iy—d—\a\ 



\D a G y {x) \ < < 



Cy^a.d I \ x \ I2 ' 

c 



c 



y.a.d 
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and this decomposition is unique. Now for any real s, the expression on the left is real. 
Therefore, taking the complex conjugate of both sides, we get 

1 «. ft, » fa f fj 

Ltft ■ L 



s m+l (! + s 2)m+l ^ ( s _ iy+i ( s + f y+i • 

Uniqueness of the decomposition then implies that for each j: 

(i) (Xj is real 

(ii) Pj = fy. 

To further characterize the coefficients, we replace s by —5 in ([3J. First, we have 

(-i) m+1 = f (-iy' +1 «y , f. (-i) J+1 A/ , f (-i) ;,+1 A/ 

and therefore 

1 _ » (-1)^ » (-1)^ " (-iy +m & 

Again using the uniqueness of the partial fraction decomposition, it follows that 

a) (-iy+-o y = o/ 

(ii) (-iy+*/3y = ft 

for each /'. The first property implies that a, = for either all odd j or all even /'. The 
second property tell us that J3y is real when j + m is even, and it is imaginary when 
j + m is odd. 

We now compute f m as the inverse Laplace transform of the sum in (|3). 

'"a- m ft '"ft- 

y=o y=o •/• y=o ■/• 

Now if m is odd, we have m = 21 + 1 and 

which can be simplified to 

AmW - E T^TTv^ 1 + E 7^^' +1 cos « - E 

7 fo I 2 ./ + 1 ) ! £0 (2; + 1 ) ! jr (2y ) ! 

Similarly, when m = 21 we have 
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Lemma 1.2. Let ^d,k be a Wendland function with d odd, and define n and m by 
d = 2n + 1 and m = n + k. Then the exact form of the Fourier transform of <t>d,k is 
found by substituting the above representations of f m into &dk{ x ) = B m fi„(r)r~ im ~ 2 . 

We can now determine properties of <tv/ # that will be required for the error bounds. 
The proof is omitted as the result follows easily from the formulas computed for f m . 

Proposition 1.3. Let d and k be non-negative integers with d odd, and considered k {x) = 
B m f m ( r ) r ^ m2 - Then (•)~ 3m ~ 2 / m is an analytic function, and 

r- 3m - 2 f m (r)=0(r- 2m - 2 ) 

as r — > oo. 

Let us now take a closer look at the 1 -dimensional case. 
Proposition 1.4. Ifk e N, then there exists Cef and h e L 1 (R) such that 

x 2k+2 ^. k {x)=B k (L + t^^ x)+C -^+h{x)j , (5) 

where h is the Fourier transform ofh. 

Proof. Since the case where k is odd is similar to the case where k is even, we will only 
prove the former. Recall that 4>j ^(x) — B k fa(r)r~ 3k ~ 2 , so let us begin by examining 
the function fa. For k = 21 + 1 and some constants aj, bj, and Cj, we have 

r-\Hr) = t ajrV + £ ^cos(r) - £ cjt**^. (6) 

7=0 7=0 ;=0 r 

We can then define an analytic function fa : R — > R by 

/*(*) = M~7*(M). 

and we will have x 2k+2 <t>i , k (x) = B k x~ k+l fa(x). Now since <J>! t G L'(R) and<f> liS .(jc) = 
B k fa(r)r~ 2lk ~ 2 , fa(x) must have a zero of order 3&+ 1 at 0, and therefore /^(jc) has a 
power series of the form 

/*(*)= £ djxK (7) 

In order to verify (0, we first need to determine some of the coefficients in (|6). 
From our previous work, we know that c/ € R, and we can find a/ and £>/ as follows. 
In the partial fraction decomposition (01, multiply both sides by s m+l (l +s 2 ) m+l . By 
substituting the values s = and s = —i, we find that a m = 1 and /3 m = (— l) m /2' n+1 , 
and therefore a, = 1 /k\ and &/ = (-l)*/(fc!2* +1 ). 

We can now finish the proof by showing that 

~ , , m ~ , ^ /I f-l) fe+1 , \ sin(x)\ 
h k (x) :=x- k+l fa(x) - f - + -_l_cosW - C/ — iij 
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is the Fourier transform of an L 1 function. Since hk{x) is identically for k = 1, we 
need only consider k > 3. This can be verified by determining that h^x) has two 
continuous derivatives in L 1 , cf. IfTTl p. 219]. Considering the representation (Q, it is 
clear that has two continuous derivatives, and the decay of these functions can be 
bounded using ([§}. □ □ 



2 General Error Estimates 

In the paper J5), the authors made use of a Green's function-type condition in order 
to determine rates of approximation, cf. [5. Property A2]. They were concerned with 
approximating functions in L p (R d ) by approximation spaces Sx ( < t ) ), where X has no 
accumulation points and hx is finite. The essence of their argument is as follows. 

Property 1. Suppose T : C^{R d ) -> C c (R d ) is a linear operator, 4> <E Lj 0C (R d ), and for 
any f &C k c (R d ) we have 

f= I T f -t)dt. (8) 

Now, given an / € C^(R d ), we can form an approximant by replacing <t> with a 
suitable kernel in dHJ. Since we are interested in approximating / by Sjf the kernel 
should have the form 

K(-,t)= £ A(t,|)*(--§). (9) 
The collection of possible kernels is restricted by requiring 

X(t) cB(t,ch x )nx 

for a constant C > and additionally requiring A (t, to be in L 1 for all |. An essential 
ingredient for deriving approximation rates is showing that <t> can be well approximated 
by K. We therefore define the error kernel 

E(x,t):=®(x-t)-K(x,t). (10) 

and assume the following property, cf. |2] Property A4]. 

Property 2. Given <t>, there exists a kernel K of the form 0) and constants l> d, k > 0, 
and C > smc/i f/iaf 

\E(x,t)\ <Ch\- d ^1 + 

forallx,t e W. 

While this approach can be used to provide estimates for some popular RBFs, e.g. 
the thin-plate splines, not all RBFs are Green's functions. Therefore our goal will be 
to extend the class of applicable RBFs. Note that in J5), the authors were concerned 
with proving error bounds that account for the local density of the data sites. As our 
goal is to extend the class of applicable RBFs, we will not discuss this more technical 



—i 
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approach; instead we will assume X is quasi-uniform, and our bound will be written in 
terms of the global density parameter hx- 

We will now show that it is possible to replace Property 1 by a condition that only 
requires <t> to be "close" to a function that satisfies this property. Specifically, we will 
consider RBFs <t> that satisfy the following. 



Property V. Let T : C* 



and G € L? 



lat- 



he a pair satisfying Property 



1, and suppose G = <$>* jJ, n + V„ where <t> <G Lj oc (M. a ), jl n is a sequence of compactly sup- 
ported, finite Borel measures with \ \pL n \ \ bounded by a constant, and V„ is a sequence 
of finite Borel measures with \ \v n \ \ converging to 0. 

The following provides an error bound for approximation by Sx where X is 
quasi-uniform, and the error is measured as 



g{f,S x m) P = inf J\fS\\ m 



for 1 < p < oo. 



Theorem 2.1. Let X be a quasi-uniform set in R . Suppose 4> is an RBF satisfying 
Property 1 ' and Property 2, and let f G C^(R d ), then 



£{f,S x {<S>))p < Chi \f\ w{LP[w i ) T) 

for any 1 < p < °°. 

Proof. First, we define a sequence of approximants to / by 



/ Tf*Hn(t)K(-,t)dt 

^ d Tf*H n (t)A(t,^dt. 



!;ex 

Note that this sum is finite due to the compact support of Tf * jx n and the conditions 
imposed on A (•,•). Now 



\f-F„\ 



< 



f Tf(t)G(--t)dt- f Tf*n n (t)K(-,t)dt 
f Tf*pn(t) (*(■- t)-K(-,t))dt + Tf *v n 

[ \Tf*Hn(t)\M--t)-K(;t)\dt + \Tf*V n 



and by Property 2, 



\f-F n \ <Ch K x d 



[jTf*Hn(t 

JR d 



h x 



-i 



dt + \Tf*v n \. 
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Therefore 



\\f-F n \\ p < Ch K x 



K—d 



Tf*pL n \*\ 1 



-/ 



\Tf*v n \ 



< c/^||r/|| p || M „|| + ||r/|| ; j|v„| 



by the generalization of Young's inequality stated in equation (HJ. Hence, for N suffi- 
ciently large (depending on X), we will have 



\\f-F N \\ p <Chl\\Tf\\ p . 
□ 



□ 



Now that we have established this result for C k c , we would like to extend it to all 
of W(L p (R ll ),T). This is accomplished in the following corollary using a density 
argument. 

Corollary 2.2. Suppose f £ W(L p (M. d ),T), then 

g(f,S x (<b))p < Chi \f\ w{LP{m i ) T] 

Proof. Let /„ € C~ be a sequence converging to / in W [LP (M. d ), T), and let F n g Sx ( < I ) ) 
be a sequence of approximants satisfying 



\\fn-F n \\ p < Ch^ \fn\ W {LP^ : 



Then 



< \\f-fn\\ p + \\fn-Fn\\ p 



< \\f-fn\ 



-Chf\fn 



W(LP(K d ),T) ■ 



Now since /„ converges to / in W(L p (M. d ), T), the result follows. 



□ 



□ 



3 Special Cases 

In this section we will provide some examples, showing how to verify the properties 
listed above. However, we will first need to take a look at local polynomial repro- 
ductions. As in 0, this will be a crucial part of verifying Property 2. To that end let 
<t> G C" (W 1 ) be the RBF under consideration, and let P denote the space of polynomials 
of degree at most n — 1 on M. d . For a finite set Y C R d , let Ay be the set of extensions 
to C(W ! ) of linear combinations of the point evaluation functionals 8 y : P\y — > R. We 
now assume the following, cf. J5] p.9]: 

(i) There is a constant C\ such that X{t) C B(t,C\h x ) HI for all t e W 1 . 

(ii) For all t , there exists A? € ^-x(t) such that Xt agrees with 8 t on P. 

(iii) 1 1 A? 1 1 < Ct for some constant C2 independent of t . 
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Based on these assumptions, X t takes the form Y.%EX{t)A(% > f )^£> an d we can define a 
kernel approximant to <I> by 



K{x,t) = X,(<P(x--)) 

= £ A(f,$)*(*-£). 

One way to verify the validity of our assumptions is the following. For each x £ 
hx^ a ', we denote by Q x the cube of side length hx centered at x. This provides a 
partition of W 1 into cubes. We then associate to each cube Q x the ball B(x,Cj,hx) for 
some constant C3, and for each t in a fixed Q XQ , we define X(t ) =X nB(xo,C?,hx). By 
choosing C3 appropriately and bounding the possible values of hx, one can show that 
there exists A, satisfying the above properties, with C2 = 2, for all t 6 Q XQ , cf. Il20l 
Chapter 3]. As xo was arbitrary, the result holds for all t e W 1 . Note that by using 
the same set X(t) for all t in a given cube Q XQ , we are able to choose the coefficients 
A (4 , f ) so that they are continuous with respect to t in the interior of the cube. To see 
this, let {pi}™ = t be a basis for P and let B(x 0l C3h x ) DX = {|y}J =1 . We now define the 
matrix M;j := Pi(£,j) and the f dependent vector j3(f); := Pi(t). Since A* reproduces 
polynomials, there exists a vector a(t) such that Ma(t) = j3(f) for all f g Q XQ . We 
could therefore choose a specific a(f) by means of a pseudo-inverse, i.e. ait) := 
M T (MM T )~ l /3(f). In this form it is clear that a is a continuous function of t and 
X t := Y!i=i a i(t)S%i satisfies properties (i) and (ii). Property (iii) follows from the fact 
that the pseudo-inverse gives the minimum norm solution. 

3.1 Sobolev Splines 

We are now in a position to prove Property 2 for the Sobolev Splines Gy where 7 is a 
positive even integer. In particular, we will show that there is a kernel Ky so that 

\Gy{x-t)-Ky(x,t)\ <Ch y x d U + ^ 

when d is odd and obtain a similar result for d even. From the argument above, we 
know that there exists Xt that reproduces polynomials of degree 7 and satisfies the re- 
quired form defined in the previous section. Using Xt, we define the kernel approximant 
to Gy by 

Ky{x,t):=Xt{Gy{x--)). 

Property 2 will only be verified for t = as all other cases work similarly. First 
suppose that ||x|| 2 > 2C\hx, then Gyis analytic in B(x,C\hx). If we let,/? be the degree 
7 Taylor polynomial of <t> at x, then we get 

\G y {x) - Ky(x,0) I = \G r (x) - R(x) - Xo(G r (x - •) -R(x- •)) I , 

and therefore 

\G Y (x)-K Y (x,Q)\<(l+C 2 )\\G Y -R\\ Ir{B{XtCihx)y 
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From Taylor's remainder theorem, we get 



\G r (x) - K Y (x,0)\ < , +l yh x WGyWwr+u^^c^)) ■ 



and applying proposition ! 1 . 1 1 gives 

\G r (x)-K r (x,0)\ < Chl +1 \x\2 d ^ (11) 
< Ch 7 f d (l + ^\ d \ (12) 

Now if | L < 2C\hx and d is odd, we define R to be the degree 7 — d — 1 Taylor 
polynomial of Gy at 0; if d is even, let R to be the degree 7 — d — 2 Taylor polynomial 
of Gy at 0. For d odd, we get 



\G r (x) - K Y (x,0)\ <Ch\ d \\G 
and by Proposition ll.il we have 

\G r (x) - K r (x,0)\ <Ch Y x d 



Y\\wY- d <°°(B(0,3C l hx)) ' 



Gy(x)-K Y (x,0)\ <Ch x + ) ■ (13) 



By making use of the assumption | |x| | 2 < 2C\hx, we finally get 

h x 

Similarly for d even, we have 

|G r (x) - K 7 (x, 0)| <Ch x - d - 1 + ^ (14) 

Putting these estimates together, equations ( fT2] i. (fT~3T >. and ( TBI imply Property 2 
for the Sobolev splines. Note that Property 1" is clearly satisfied for G — Gy and T = 

(l-A)r/ 2 . 

Theorem 3.1. For any positive even integer 7 > d + 1 and any 1 < p < °° the following 
rates are valid for approximation by the Sobolev splines Gy 

i) If d is odd 

&(f,Sx(Gy)) p < Chi \f\ Lrp{R d ] 

ii) If d is even 

<£{f,Sx(Gy)) p < Ch y x X |/| £ r,P(E rf ) ■ 
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3.2 Wendland Functions 

We now come to the reason for studying perturbations: to derive approximation results 
for the Wendland functions. We will begin by using the polynomial reproducing func- 
tionals X t to verify Property 2. Note that <J> ( /.jt(jc — ?) and Xt(<&d,k(x — •)) me b° tn zero 
for | |jc — 1\\ 2 > 1 +C\hx- Therefore in order to verify Property 2, it suffices to show 
that for \\x — t\\ 2 < 1 + C\h x we have \<Sf^{x — t) — hz{$>d,k{x— 0)1 ^ Ch^. Similar 
to the Sobolev splines, we will only verify this inequality for t = 0. To begin, fix x, and 
let R be the 2k— 1 degree Taylor polynomial of <t> f / j- at x. Then 



|«& d) *(x)-Ao(4> d) *(x-0)| = |(®^-*)W-^)((*^-*)(*-0)| 



< ||Ao||||*rf,t-/?| 



We have therefore shown that each <£></ 1 satisfies Property 2 with fc = 2^:. 

The challenge in proving rates for the Wendland functions appears when trying to 
prove Property 1'. To accomplish this, we can work in the Fourier domain, where 
the convolution becomes a standard product. For example, to show that G = <i>* jX for 
some jj, € M(R d ), we can verify that G/<i> is the Fourier transform of some fi 6 M(W ! ). 
The difficulty lies in characterizing the space of Fourier transforms of M(W ! ). This is 
known to be a very difficult problem, cf. 0. However, in certain situations, we are able 
to make this determination. Our approach will be to work under the assumption that the 
functions being approximated have added smoothness, and we will additionally work 
in odd space dimension d, since we know the exact form of the Fourier transform of 
<J>d £. In this situation, we can choose our Green's function G to be the Sobolev spline 

G2d+2k+2- 

Lemma 3.2. Let G and <i>d,k be defined as above with k € Z + . In this case we have 

IX := (G/4>) v GL 1 ^) 



Proof. First, it is clear that G/4> is a radial L 1 (W 1 ) function, so Li e L°°(]R rf ). Addition- 
ally, we can write 



1" 



where J is a Bessel function of the first kind. In this form, we can bound jj. for r large, 
and hence show that it has sufficient decay to be in L l (M. ). Verifying the decay of jj, 
is accomplished using the fact that 

^-{z v J v (z)}=z v J v -i(z) 
dz 

and integrating by parts (d + 5)/2 times. Note that Proposition |1.3| is used to bound the 
decay at each step. □ □ 
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With pi defined as in the previous lemma, we get 

Gid+ik+i = <&d,k * M > 

and by choosing \i n to be the restriction of fi to B(0,n), we can see that Property 1" is 
satisfied. Therefore we have the following theorem. 

Theorem 3.3. For any positive integer k and any odd space dimension d, we have 
£(f,Sx(®d,k))p<Chf\f\ L u + u + 2, P 

for 1 < p < oo. 

Upon inspection, one can can see that we have to assume a higher order of smooth- 
ness on the functions being approximated than the rate of approximation. Notice 
that the approximated functions / are assumed to be in ]_, 2d + 2k + 2 'P ) but the rate of 
approximation is only h 2 ^. Despite this shortcoming, the derived rate is better than 
what is currently known in certain cases. In lEUl . Wendland determines a rate of 
^(r/+2£+i)/2 d(i/2 l/p)+ f or bounding the L p (£l) error in approximating functions from 
W( d+2k + l )l 2 - 2 (a), where Q. is a bounded domain in W l . First, notice that when approx- 
imating arbitrarily smooth functions and 1 < p < 2, Theorem l3 . 3l provides a better rate 
as long as k > (d+l)/2. However, the improvement is most noticeable when p is close 
to oo. Also note that the prior rate suffers the same drawback as Theorem [33j when 
p = 1, the L l smoothness needed to be in W^ d+2k+1 ^ 2 ' 2 (Q.) is greater than d + 2k+ 1. 

3.3 Improved Estimates for the Wendland Functions 

In the generality of Theorem 13. 3 1 the assumed smoothness is greater than the rate of 
approximation. One source of this disparity was that when verifying Property 1', we 
compared <P c i,k to the Green's function God+ik+i which has a Fourier transform with a 
faster rate of decay. In this section, we propose an alternate method and show that for 
d = 1 we have <t>i ^ = Gik+2 * M f° r some finite Borel measure ji and hence improve 
the error estimates. Notice that the Fourier transforms of <t>i ^ and Gik+2 have the same 
rate of decay, so in this situation, we are no longer losing orders of approximation from 
the perturbation technique. 

Our approach to verifying Property 1' will be to view M(R) as a Banach algebra 
and use an invertibility result of Benedetto. The first step is to show that <t>/G is the 
Fourier transform of some jj. E M(R). Then the following theorem will give condi- 
tions for /I -1 = G/O being the Fourier transform of an element of M(R). Recall the 
decomposition of measures: a Borel measure ji can be written as \i a + fi s + fij. Using 
this notation, we state the following theorem of Benedetto. Note that the author proved 
this theorem in a more general setting with M replaced by an arbitrary locally compact 
Abelian group. 

Theorem 3.4. ([3, Theorem 2.4.4]) Let ji e M(R) such that \fl \ never vanishes and 

\\HsW < inf \fid(x)\. 

xeM 

Then £i~ l is the Fourier transform of an element o/M(M). 
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We now fix k > 1, and we must select a Green's function G and corresponding 
differential operator T so that <&\ j,/G satisfies the necessary conditions. Considering 
the decay of 4> u , we choose G= (1 + \-\ 2k+2 )- 1 and T = (l+ (-l)* +1 A fc+1 ). Note 
that for 1 < p < oo we have W(L P (R),T) is equivalent to the spaces W 2k+2 ' P (R) and 
L 2k+2 >P(R). 

With this Green's function G, we can verify the hypotheses of Theorem 13.41 By 
Proposition |1.4| we know that for some > and C £ M. we have 

$1*^ - /I f-l)* +1 , , sin(x) f/ ,\ 

-^(x) = * lik+ 2?» + cos(x) - C-±±+h(x)) 

where h is the Fourier transform of an L 1 function, h. The fact that <t>i ^/G is positive 
follows from the positivity of 4>i,t, cf. [20, Chapter 10]. Now let fi be the measure 
defined by 

M = + v^B* (I* + ^^(5-, + 50 - ^ hu] + ^Lfc) - 

Then is a finite Borel measure with /i = <t>i ^/G. Additionally, has no singular 
continuous part, and 

Therefore Theorem 13. 4l implies that \ / fx ~ G/<$>i.k is the Fourier transform of some 
finite Borel measure p.. By letting p. n be the restriction of p, to B(Q,n), we can see that 

G = 4> u .*£ 
and hence corollary 12. 2l applies. 

Theorem 3.5. Let <t>i ^ fee as above, and let X be a quasi-uniform subset o/R. Then 
forf € L 2A:+2 <' 5 (R) and\<p<™, we have 

£{f,s x {<b hk )) p <ch 2k \fy k+2 ,,, m 

By examining this example, we can determine what the analogous result would be 
for the remaining Wendland functions. If k > 1, then it is known that 

ci(l + \\co\\ 2 2 )-^ +2k+ ^l 2 < $ d>k (a>) < c 2 (l + \\o)\\lr {d+2k+1)/2 

for some positive constants c\ and c%, cf. ll20l Theorem 10.35]. Therefore, we could 
choose G to be the Sobolev spline of order d + 2k + 1, so that G and 9^ have similar 
decay. In this situation, G/<& ( i,k is a continuous function that is bounded above and 
bounded away from 0. If we could show that G/^ c i,k is the Fourier transform of an 
element of M(R d ), then we could apply corollary 12. 2 I to obtain the following. 

Conjecture 3.6. Let <t>d.k be os above with k > 1, and let X be a quasi-uniform subset 
ofW 1 . Then for / € L cl+2k+l -P{W') and\<p<™, we have 

S{f-,Sx{^d.k))p < Chf \f\ L d+u+i,p(jgd) ■ 
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